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We study some mesoscopic properties of electron transport by employing one-dimensional chains 
and Anderson tight-binding model. Principal attention is paid to the resistance of finite-length chains 
with disordered white-noise potential. We develop a new version of the transfer matrix approach 
based on the equivalency of a discrete Schrodinger equation and a two-dimensional Hamiltonian 
map describing a parametric kicked oscillator. In the two limiting cases of ballistic and localized 
regime we demonstrate how analytical results for the mean resistance and its second moment can 
be derived directly from the averaging over classical trajectories of the Hamiltonian map. We also 
discuss the implication of the single-parameter scaling hypothesis to the resistance. 



I. INTRODUCTION. 

In recent years, the study of one-dimensional (ID) 
tight-binding models with diagonal disorder has led to 
new important results. The growing interest in the ID 
disordered models is mainly due to two reasons. First, it 
was found that specific correlations in (random) poten- 
tials give rise to delocalized states. Early demonstrations 
of this effect are related to the so-called dimer model [1] 
that is specified by short-range correlations. It was shown 
that for particular discrete values of the electron energy 
the localization length diverges, in contrast to a common 
belief that any randomness in the potential leads in ID 
geometry to an exponential localization of all eigenstates. 
Recently, delocalized states have been experimentally ob- 
served in random dimer superlattice [2] . Similar effects of 
delocalization due to correlated disorder have been also 
predicted for some special ID models, see [3] and refer- 
ences in [4]. 

Further study of ID random models with correlated 
potentials produced even more exciting results. It was 
shown [5,6] that specific long-range correlations may re- 
sult in an appearance of the mobility edges. This means 
that a continuous range of energy arises on one side of 
the mobility edge, where the eigenstates are extended. In 
Ref. [6] the method for constructing such correlated ran- 
dom potentials has been proposed. Using this method, 
one can relatively easily construct specific random po- 
tentials that result in energy bands of a complete trans- 
parency for finite samples. The position and the width 
of such windows of transparency can be controlled by the 
form of the binary correlator of a weak random potential. 

The role of long-range correlations has been studied 
in details for the tight-binding Anderson-type model [6] , 



and for the Kronig-Penney model with randomly dis- 
tributed amplitudes [7] and positions of delta-peaks [4]. 
These results have been also extended to a single-mode 
waveguide with random surface profiles [8]. The pre- 
dictions of the theory [6] have been verified experimen- 
tally [9] , when studying transport properties of a single- 
mode electromagnetic waveguide with point-like scatter- 
ers. The latter were intentionally inserted into the waveg- 
uide in a way to provide a random potential with slowly 
decaying binary correlator. Very recently [10] the ex- 
istence of mobility edges was predicted for waveguides 
with a finite number of propagating channels (quasi- ID 
system) with long-range correlations in a random surface 
scattering potential. 

Another reason for the revival of the interest in ID dis- 
ordered models is due to the revision of the famous sin- 
gle parameter scaling hypothesis (SPS) [11] for transport 
characteristics of disordered conductors. As is shown in 
[12], the assumption of a random character of fluctua- 
tions of the finite-length Lyapunov exponent, originally 
used [13] in order to justify the SPS, turns out to be 
incorrect. Specifically, it was found that in the vicinity 
of the band edges the SPS is violated [12]. This result 
is important both from the theoretical and experimental 
viewpoint (see discussion and references in [12]). 

One of the effective tools to analyze ID tight-binding 
models is the one based on the Hamiltonian map (HM) 
approach that has been developed in Refs. [14-16]. The 
key point of this approach is a transformation that re- 
duces a discrete ID Schrodinger equation to the classi- 
cal two-dimensional Hamiltonian map. The properties 
of trajectories of this map are related to transport prop- 
erties of a quantum model. The geometrical aspects of 
the HM approach turn out to be helpful in qualitative 
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analysis, as well as in deriving analytical formulas. 

With this method many known results can be easily 
obtained for the localization length l^. In particular, 
a rigorous derivation of was made in a general form 
for energy close to the energy bands of the standard An- 
derson model with uncorrelated diagonal disorder [16]. 
This method has been also used in [6] to derive a gen- 
eral expression for the localization length for any weak 
potential. All these studies refer to infinite samples and 
are based on the evaluation of the infinite-length Lya- 
punov exponent that is determined by the evolution of 
the classical Hamiltonian two-dimensional map. 

In this paper we present further developments of the 
HM approach, that allows us to obtain important trans- 
port characteristics in ID random models. The main at- 
tention is paid to the resistance of finite samples that are 
described by the standard tight-binding Anderson model 
with weak white-noise potential. We derive analytical 
expressions for the mean value of the resistance and its 
second moment. These are the quantities that are easy 
to measure experimentally. We also discuss the implica- 
tion of the SPS hypothesis for the resistance in strongly 
localized regime. 

II. HAMILTONIAN MAP 

The stationary discrete Schrodingcr equation for ID 
tight-binding models with diagonal disorder can be writ- 
ten as follows, 

V>„+1 + tpn-1 = (E + € n )%p n , (1) 

where E is the energy of an eigenstate and e„ is the site 
potential. It is known that this equation is equivalent to 
diagonalization of a tridiagonal matrix [17] or to a clas- 
sical Hamiltonian map for canonical variables x n = ip n 
andp„ = (ip n cos fj,-ip n -i)/ sin fj, (sec, e.g. [14]), 

( x n+ i \ = ( (cos/i + A n sin/i) sin/z \ / x n \ 
\Pn+i J V (Ai cos /z- sin/z) cos^i J \p n J 

The canonical variables x n and p n can be considered as 
position and momentum of a linear oscillator subjected 
to linear periodic delta kicks. The amplitude of the nth 
kick depends on the electron energy at the nth site poten- 
tial, A n = 6 n / sin/z. In this approach, the amplitude 
ip n of an eigenstate at the nth site is given by the posi- 
tion of the oscillator at time t n = n. Therefore, global 
properties of the eigenstates can be studied by exploring 
the time dependence of the classical map (2). The energy 
E is related to the angle /z via E = 2 cos /z. 

One can see that the representation (2) is, in essence, 
the Hamiltonian version of the standard transfer matrix 
method. Indeed, starting from two initial values ipo and 
ip-i one can compute ip n and ip n -i according to the re- 
cursion given by Eq.(l) or Eq. (2). Such a representation 



allows one to determine the Lyapunov exponent A by 
running the trajectory of the map Eq. (2) from n = 
to n — > oo. According to the standard definition (see, 
e.g. [18]), the inverse of the Lyapunov exponent A _1 (_E) 
gives the localization length (E) of the eigenstate. 

In many aspects the Hamiltonian representation (2) 
is more convenient than the standard one based on the 
original equation (1). The effectiveness of the treatment 
of the map (2), instead of (1), is clearly manifested if the 
action-angle variables (r, 6) are introduced according to 
the following transformation, x — rsin# and p — rcos9. 
Then the map Eq. (2) takes the following form, 

sin 6 n+1 = D' 1 [sin (6 n - /z) - A n sin 6 n sin /z] , 

(3) 

cos 9 n+ i = D n 1 [cos (6 n - /z) + A n sin 9 n cos /z] 
where 

D n = — = Jl + A n sin (29 n ) + Alsin 2 9 n . (4) 

Using action-angle variables (3), (4), the inverse lo- 
calization length is written in a quite simple form (see 
details in [16]), 

J- 1 = A= lim IV>|%±i|= lim IVln^i 

°° L^oo L ^ ip n L — 00 L r n 

71=1 Tn 71=1 

= i (In (1 + A n sin 26 n + A 2 n sin 2 6 n )) n • (5) 

Here the brackets (•••)„ stand for averaging over n, i.e. 
along the trajectory of the map (3), (4). It is important 
to note that the above expression for A depends only on 
the angle 6 n and not on the radius r„. The above rela- 
tion is correct for the energies E not very close to the 
band edges E = ±2 where /i = 0, tt. At the band edges, 
there is additional contribution to the localization length 
that depends on the ratio sin# n+ i/sin#„ [16]. In this 
paper we consider the case when the energy E is inside 
the allowed band, \E\ < 2. Moreover, we also exclude 
the band center E = where the localization length has 
a singular behavior that requires specific treatment (see 
[16] and references therein). 

Apart from the abovementioned restrictions, the rela- 
tion (5) is valid for any potential e n . However, analyti- 
cal treatment of the localization length is possible in the 
two limit cases of a weak or strong potential. Since in 
the following we consider the case of weak disorder only, 
| A n |<C 1, let us demonstrate how the localization length 
can be derived from (5). We specify that the distribution 
function of site energies e„ is flat, P (e„) = 1/W, within 
the region | e n |< W/2, with the variance (e 2 ) = W 2 /12. 
Then, for W <C 1 the logarithm in (5) can be expanded, 
and by taking into account that in the lowest approxima- 
tion with disorder the values of A n and 9 n are statistically 
independent, one can easily obtain, 
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For the first time this expression was derived by Thou- 
less [19]. It works quite well over the whole range of 
energies apart from the band center E = 0. The analyti- 
cal treatment of the expression (5) was the main interest 
for different cases, including correlated disorder. Unlike 
our previous studies [4,6,7,14,16], in what follows we ad- 
dress a question about the resistance of finite samples of 
size L. Therefore, we are interested in the properties of 
the map (2) on the finite time scale, for n = 1, L. 

Our further consideration is based on the expression 
for the transmission coefficient Tl in terms of classical 
trajectories of the map (2) [15], 
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Here t\,l and r 2j L stand for the radii of the two compli- 
mentary trajectories obtained by iterating the map (3) 
up to the last site of the sample, n = L. Each of the tra- 
jectories is specified by the initial conditions at L = 0, 
namely, n,o = T2,o = 1) #i,o = 0, and 6 2 ,o = n/2. 

As one can see, all statistical properties of the trans- 
mission are entirely determined by the evolution of the 
two complementary trajectories of the classical map. 
Some analytical and numerical analysis of the expression 
(7) have been recently performed in Ref . [20] . The ques- 
tion under study was the statistical distribution of the 
transmission coefficient Tl- In particular, it was shown 
that the correlations between the two classical trajecto- 
ries are different for the ballistic and localized regimes, 
giving rise to the different distribution functions of Tl ■ In 
next sections we address the question of global properties 
of the resistance Rl — T£ , by paying main attention to 
the mean values (Rl) and (R L )- 



III. RESISTANCE. 

A. First moment. 

According to Eq. (7) the resistance Rl = T^ 1 of a 
sample of length L is given by the following formula, 
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Here the final radii rf L (i = 1,2) are expressed through 

the kicks amplitudes A n and the phases On = 0i tTl of the 
two trajectories. Using Eq. (4) we get, 

L-l 

r h = LI i 1 + A n Sin ' 9 n + A " Sin W n) ' ( 9 ) 

n=0 



It is convenient to represent the mean value of rf L in the 
equivalent form 



'L-l 



(rf ;i ) = / cxp ]T In (l + Al sin 2 6»« + A n sin 29®] 



„ n=0 



Then, in the limit of weak disorder, | A n |<C 1, one can 
expand the logarithm in Eq. (10). Separating quadratic 
and linear terms, we represent (rf L ) in a form of the 
product of two factors, 



'L-l 



L-l 



(rU) =(exp E Al sin 2 «#> _ I £ A 2 sin 2 2^ 



, n=0 
f L-l 



n=0 



, n=0 



The first factor contains the sums X)n=o ^-n sm2 
where (j$ stands for either 6$ or 2#. Since L » 1 
these sums are self-averaged quantities, with small vari- 
ance. Thus, they can be substituted by their mean val- 
ues, 

E A ' ^ ^ W ) = 4 A , (12) 
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n=0 
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Therefore, the first factor of Eq.(ll) takes the form, 



'L-l 



L-l 



exp E A l Sin2 W + E A n Sin2 W n M « <*P (2A) 



.71=0 



n=0 



(14) 



The second factor in Eq.(ll) exhibits strong fluctua- 
tions. We can calculate its mean value using the distri- 
bution function V(SX), 



/oo 
f(S A )V(S A ) dS A 
-oo 



(15) 



where /(Sa) = exp(SA) with Sa = E^o An sin 26"^ l \ 
The variable Sa is a sum of random independent num- 
bers, therefore, according to the central limit theorem 
P(Sa) is a Gaussian function with variance 4A. The dis- 
tribution of Sa is shown in Fig. 1 and one can see quite 
good fitting of the numerical data by a Gaussian func- 
tion. As a result, the second factor of Eq.(ll) can be 
calculated explicitly, 

1 f°° 

(exp (S A )) = -= / exp (S A - SA^Si) dS A 

Vo7rA J-oo 

= exp (2A) . (16) 
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FIG. 1. Histogram for the probability distribution of 
Sa = J2tZl A ^ sin2e n } plotted for E = 1.5, W = 0.1, 
L = 4200 and A = L/l^ = 1 for 10,000 realizations of the 
random potential e„. The smooth curve is a Gaussian distri- 
bution with variance a 2 — 4. 



d2 _ 1 . 1 2 , 1 2 , 1 4 i 4 , 1 2 2 

L ~ 4 4 ri ' L 4 r2 ' L 16 ri,L Hf 2 ' 1 - 8 ri ' i? ' 2 ' i ■ 

(19) 

One can see that apart from the second and fourth mo- 
ments of ri t L, this expression contains a product r\ L r\ L . 
The mean value of this term depends on the correlations 
between the two complementary trajectories of the clas- 
sical map (3). This fact strongly complicates analytical 
treatment. In what follows, wc consider the two limiting 
cases of the ballistic and localized regimes separately. 

First, let us start with the evaluation of the mean value 
of the fourth moments rf L . This can be done by using 
Eq.(9) in the same way as described above, directly and 
without restrictions, 

(ft) - ( ex p f 2 E A « sin2 W - E A * sin2 2e ") ) 

\ V n=0 n=0 / / 

x ^exp ^2 A n 20<?^ ^ . (20) 



Substituting Eqs.(14) and (16) into Eq.(ll), we obtain 
the following expression for the mean value of rf L , 

<r 2 i )=cxp(4A) . (17) 

This leads to the final formula for the mean value of the 
resistance (8), 

(i? L ) = i[l + exp(4A)] . (18) 

Note that this expression is valid for any value of the con- 
trol parameter A. It gives the resistance in the metallic 
(A <C 1) and insulator regimes (A 3> 1) as well as the 
crossover (A ~ 1). The only condition that we have used 
here is the condition of weak disorder. 

In the metallic regime, the fluctuations of the transmis- 
sion coefficient are weak and it is a self-averaged function. 
Therefore, the mean value of conductance (G) is the same 
as 1/ (R). However, in the localized regime, both the con- 
ductance and resistance are not self-averaged functions 
that lead to distinct values of (G) and 1/ (R). Indeed, 
it is well-known that the mean conductance has a pre- 
exponential factor <~ (L/^oo) 3 / 2 . Derivation of this factor 
requires some special efforts [18]. This can be under- 
stood from Eq. (7) which contains the random radii in 
the denominator. Unlike this, in the case of resistance 
(8) these random variables reside in the nominator and 
strongly simplify the calculations. 

B. Second moment. 

From Eq.(8) the second moment of the resistance can 
be expressed as follows, 



Here we have used the weak disorder condition A^C 1 
and have kept only terms up to O (A 2 ) . We have also 
presented the mean value in the form of the product of 
two factors, one with a small variance (first factor) and 
another that reveals strong fluctuations (second factor). 
Inside the first factor, the sum X)«=o A n sm2 with 
as 9n^ or 26n\ is a self-averaged quantity. There- 
fore, after substitution by its mean value we get, 

/exp U J2 A\ sin 2 0« - ]T ^l 2 sin 2 20<A ) « cxp (4A) . 

\ \ n=0 n=0 / / 

(21) 

For the second factor of Eq. (20) we proceed in the same 
way as we did in the derivation of Eq. (16). We use 

the fact that the quantity Sa — X)«=o ^™ s ^ n ^™ 1S dis- 
tributed according to the Gaussian (see demonstration in 
Fig.l) with the variance 4A. Therefore, the mean value 
of the second factor in Eq. (20) is given by, 

1 I" 00 

(exp (2S A )) = -= / exp (2S A - SA^S*!) dS A 

v87rA J-oo 

= cxp(8A) . (22) 

Substituting Eqs. (21) and (22) into Eq. (20) we obtain 
the following result for the mean value of rj L , 

(rf >L ) =exp(12A) . (23) 

Let us now consider the term r\ L r\ L . We can write 
the mean value of this quantity in the following exact 
form, 
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^exp |^ln(l + ^ sin 2 0™ + A n sin 26' r ( 1 1) ) 
x exp J ^ In (l + A 2 n sin 2 ^ 2) + A n sin 2<9, (2) ) 



(24) 



. n=0 



Using the weak disorder condition we can expand the log- 
arithms and present the correlator (r\ as follows, 



Ml) = ^ cx p|E A ' Z «| 



'L-l 



X ( CX P 1 E An ( S[n2d n ) + S ' m20 n ) ) 



. n=0 



«exp{L(,4 2 )(Z„)} x (exp(S P )) 
Z n = (sin 2 0W + sin 2 ^ 2) ) - I (sin 2 20« + si: 



(25) 
(26) 



where 



2 2 0(2) 



Here we introduced a random variable Sp, 

L-l 

s P = 



^2 A n (sin 



29^ 



sm 



n=0 



(27) 



(28) 



and splited the mean value (r\ L r\ L ) into two factors, 
one with small variance and another with strong fluc- 
tuations respectively, see Eq.(25). Now, we have to take 
into account the correlations between the phases 9^ and 

(2) 

9„ . In the first factor we have the following sclf-averaged 
quantities: ^n=0 A n-> E«=o ( sin2 ^ + sin2 ' and 



J2nZo ( sin2 20 n ] + sin 2 2<4 2) ) . The mean value of (A 2 ) 



g(2) 



is 8/Zoo, while the mean values ^sin 2 ^ ; +sin 2 9 { n>j and 

( sm sin 2 29 n^) are very close to 1. Substituting 

these mean values in the first factor of Eq.(26) we get, 

(A. L r\ L ) = cxp (4A) (cxp (S P )) . (29) 

Now, we have to take into account the correlations be- 
tween the phases 9^ and 9^ . These correlations is the 
main problem for further analytical treatment. It turns 
out that the correlations between phases 9„^ and 9 n 2 ^ 
are very different in the ballistic (A <C 1) and localized 
(A>1) regimes. 



,2/3(2) 



1. Ballistic regime. 

In the ballistic regime the exponential cxp Sp is close 
to one, and we expand it up to quadratic terms, 



<exp(S P )> « I + 5 2 



E^sin^W 

\n=0 



+ 



^ Al sin 2 2tf! 2) 



Here the term Si has the following form, 
S, = ( >> 2 sin (2^) sin (20^) 




(30) 



(31) 



This term describes the correlations between the phases 
9^ and 9n ^ of the two classical trajectories that start 
from complementary initial conditions. It should be 
stressed that the correlation term keeps track of the 
change in the correlations along a finite sample, when 
running from n = to n = L — 1. Taking into account 
that for the white-noise sequence e n the fluctuations of 
A n and 9 n are statistically independent [6], we can ex- 
press Si as follows 



L-l 



S 2 = 8l^J2MK). 



n=0 



Here we introduced a binary correlator i? 2 
R 2 - n ( 20\; ) sin ( 20?)' 



(32) 



(33) 



that is expected to depend on the scaling parameter 
A„ = n/loo only. Note that the average in Eq.(32) is 
taken over different realizations of the disordered poten- 
tial for a fixed value of n = 1,...,L — 1. Due to the 
ergodicity, this average is performed over different real- 
izations of e„ after running along the trajectory up to a 
fixed n. 
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FIG. 2. Numerical data for the correlator R2 versus scal- 
ing parameter A„ = n/lac is shown for the transition from the 
ballistic to localized regime for E = 1.5 and W = 0.1. The 
average was done over 10 4 realizations of the disorder with an 
additional "window moving" average, in order to reduce large 
fluctuations. 
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A specific character of the correlator R 2 is clearly seen 
from the data reported in Fig. 2. As one can see, the 
correlator R 2 changes from —1/2 in the ballistic regime, 
X n -C 1, to 1/2 in the localized regime, A n > 1. Tak- 
ing into account that (sin 2 n ) = 1/2, the latter means 
synchronization of the phases in the localized regime, 
6n^ 9n \ sec details in [20]. In the ballistic regime the 
phases are anti-synchronized, On ~ —On \ however, this 
holds only for very small values of A as one can see from 
the graph in Fig. 2. This behavior is a manifestation of 
peculiar statistical correlations between the phases On" 1 
and On- 

In order to evaluate analytically the correlator R 2 in 
the ballistic regime, we use the approximate map for the 
phase n that can be directly obtained from Eq. (3) in 
the limit of weak disorder, 

a a , sin 2 n 

n +i =0 n -fi + e n — . (34) 

sin 

Iterating this map, one can express an angle n in terms 
of the amplitudes eg, e±, e n -i of all previous kicks for 
a fixed value of \i (energy E). Taking into account the 
initial conditions 0^ = and 0^ = n/2, the following 
explicit formulas are derived for the phases, 

n X) ~ -nfj, + sin -1 fi J2i=o e t sin 2 (z^) 

- sin~ 2 nJ27=o 6i sm3 (*M) cos(i^) (35) 

- sin~ 2 fi Yn<k e i e k sin 2 (i^) sin(2fc^) + O (e 3 ) , 
and 

On^ « § - n^i + sin -1 fJ,J2i=o e « c °s 2 (*m) 

- sin~ 2 fi Y^iZa e i cos3 (i(J>) sin(i/x) (36) 

- sin -2 nY%<k €i£k cos2 (*M) sin(2fc/i) + O (e 3 ) 

We substitute the above expansions into Eq. (33) for 
the correlator R 2 . Neglecting the contributions of the 
fast oscillating terms sin m (kn/j,) with positive integers k 
and m, after some algebra we get, 

Ri{X n ) = -\+ 4A„ - 16A 2 + O (A 3 ) . (37) 

Numerical data in Fig. (2) confirm the validity of this 
estimate in the ballistic regime, A n < 1 (see also [20]). 
Then the two-point correlator S2 takes the form, 

^ = -4A + 16A 2 + O (A 3 ) , (38) 

Now we come back to Eqs. (29), (30), and (31). The 
mean values of the two sums in Eq. (30) can be easily 
evaluated and for the correlator (r\ L r\ L ) we obtain, 

{tIltIl) = [1 + 16A 2 + O (A 3 )] cxp (4A) . (39) 



Now the second moment of the resistance in the bal- 
listic regime can be expressed as, 

(R 2 L ) = 1 [2 + (5 + 16A 2 ) cxp (4A) + cxp (12A)] . (40) 

Since this expression is valid in the quadratic approxima- 
tion with A 2 , we expand the exponential terms and obtain 
the final formula for the second moment of resistance in 
the ballistic regime, 

(R 2 L ) = 1 + 4A + 16A 2 + O (A 3 ) . (41) 

There is a good agreement between this analytical result 
and numerical curve shown in Fig. 3 up to A w 0.3. 




1 , 1 1 1 , 1 

0.0 0.1 0.2 0.3 

X 

FIG. 3. Plot of Eq. (41) (solid line) against numerical data 
(open circles) within the ballistic regime. Numerical data 
are obtained for E = 1.5, W = 0.1, l x = 4200 and differ- 
ent lengths L, with the average over 50, 000 realizations of 
disorder. 



2. Strongly localized regime. 

Let us now consider strongly localized regime, A > 1. 
As was mentioned, in this regime the phases fluctuate 
coherently, i.e. o'iP ~ 0^ ■ This property helps to evalu- 
ate the mean value (Sp), see Eq. (28). Using the same 
procedure that was used for the ballistic regime, we find 
that in the localized regime the distribution function of 
random variable Sp is the Gaussian with zero mean value 
and variance 

(S 2 P ) - (Sp) 2 - 16A . (42) 
Numerical data in Fig. 4 confirm this analytical result. 
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FIG. 4. Numerical data (solid line) for the distribution of 
Sp plotted against the Gaussian (dotted line) with the zero 
mean and variance a 2 = 16. The histogram for the distri- 
bution of Sp fits well by a Gaussian curve. Numerical data 
are shown for E = 1.5, W = 0.1, L = 42000 and L/l x = 10, 
with additional average over 10 4 realizations of the random 
potential. 

The mean value of the normally distributed random 
variable Sp can be easily evaluated, 

(exp (Sp)} = cxp (8A) . (43) 

Then, substituting Eq.(43) into Eq.(29), we get the an- 
alytical expression for the mean value of the correlator 

(jlAx) = (12A) . (44) 

In the localized regime this correlator is exponentially 
large. This means that the trajectory of the map Eq. (2) 
extends far away from the origin [14]. Now it is easy to 
show that the variance R\ also grows exponentially in 
the localized regime, 

( J R!)^iexp(12A) + icxp(4A) + i. (45) 

Direct numerical computation of the mean value (R L ) 
gives rise to serious problems, since in a strongly local- 
ized regime the variance as well as all the moments of R\ 
fluctuate enormously with realization of the disordered 
potential. This is a reflection of the well known fact that 
Rl (as well as the conductance Tl) is not a self-averaged 
quantity. Instead, it is more physical to make an average 
of the logarithm of the moments. 

IV. DISCUSSION AND CONCLUDING 
REMARKS 

To the best of our knowledge, there is no systematic 
study of the statistical properties of the resistance Rl 



for the tight-binding models. Only brief discussion of the 
moments of Rl can be found in Ref. [21], where the anal- 
ysis was done with the use of the "two-scale approach". 

In the spirit of the discussion of the single parameter 
scaling [12] of the conductance, it is also interesting to 
consider the first and second moments of ln(i?g). Due 
to the simple relation Rl = T^ 1 between the resistance 
and conductance, it is clear that in the strongly localized 
regime, where the SPS is assumed to work well, the re- 
sults [20] obtained for Tl are also valid for Rl- Indeed, 
since ln(i?^) = — ln(Tt), the following relations hold 

(lni? L ) = 2A; (in 2 R L ) = 4A + 4A 2 , (46) 

Var (lnR L ) = (in 2 R L ) - (lnR L ) 2 = 2 (lni? L ) . (47) 

The last formula corresponds to the single parameter 
scaling for resistance. We remind that the SPS is ap- 
plied for conductance Tl [12] and results in the relation, 

Var (lnTj,) = -2(lnT L ) . (48) 

In conclusion, in this paper we have analytically stud- 
ied the main properties of the resistance Rl of finite ID 
samples with random potentials. Our analysis is based on 
the Hamiltonian version of the transfer matrix method. 
The main equation of this method is the Hamiltonian 
map that replaces the discrete Schrodinger equation. In 
this way all statistical properties of resistance can be ex- 
pressed through two complimentary trajectories of this 
map. Statistical correlations between these trajectories 
play an essential role and an analytical evaluation of the 
moments of the resistance is possible in the situations 
when the binary correlator of phases is known. At the 
same time, the mean value of Rl is independent of this 
correlator. Therefore, the expression for (Rl) can be 
found exactly, in the whole range of the scaling param- 
eter, < A = L/loo < oo. In contrast, the mean value 
of the second moment, (-R 2 ), depends on the binary cor- 
relator i?2- Analytical evaluation of this term was found 
to be possible in the limiting two cases of ballistic and 
localized regimes. The derived analytical expressions for 
the resistance and its second moment correspond quite 
well to the numerical data obtained in a wide range of 
the model parameters. 
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